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Using tools of nonequilibirum mechanics, we study a model of self-propelled hard rods on a 
substrate in two dimensions to quantify the interplay of self-propulsion and excluded-volume effects. 
We derive of a Smoluchowski equation for the configurational probability density of self-propelled 
rods that contains several modifications as compared to the familiar Smoluchowski equation for 
thermal rods. As a side-product of our work, we also present a purely dynamical derivation of the 
Onsager form of the mean field excluded volume interaction among thermal hard rods. 



I. INTRODUCTION 



Self-propelled particles draw energy from internal or external sources and dissipate this energy by moving through 
the medium they inhabit. A wide class of systems, including fish schools, bacterial colonies, and monolayers of 
vibrated granular rods can be described within this paradigm. These systems exhibit rich collective behavior, such 
as nonequilibrium phase transitions between disordered and ordered (possibly moving) states and novel long-range 
correlations and have been the subject of extensive theoretical numerical 0-01 and experimental investigations 

[vl-fioj in recent years. 

Self-propelled particles are elongated in shape and have a self-replenishing momentum along one direction of their 
long body axis. They generally experience attractive and repulsive interactions, both of a direct nature and mediated 
by the medium. One generic interaction that is relevant to all self-propelled systems is the short-range repulsive 
interaction arising from the finite size of the self-propelled units. Our goal here is to understand the interplay of 
self-propulsion and short range repulsive interactions in controlling the collective dynamics of the system. To this 
end, we consider the simplest implementation of short range steric repulsion, which is the hard particle limit, and 
consider a minimal model of self-propelled hard rods moving on a substrate in two dimensions. We will show that 
self-propulsion modifies the momentum exchanged by hard rods upon collision and the resulting mean-field excluded 
volume interaction as compared to the Onsager result for passive rods [H, [l^ • 

The simplest model of equilibrium nematic liquid crystals is a collection of long, thin hard rods In the 

overdamped regime and at low density, the dynamical properties of the system are described by a Smoluchowski 
equation for the configurational probability density, c(ri,Ui,t), of finding a rod with center of mass at ri and long 
axis oriented along the unit vector Ui at time t, given by 



+ Vi • J + 7^l • 

at 







(1) 



where TZi = Ui x dt^^ is a rotation operator and J and J'^ are translational and rotational currents, respectively, given 
by 



Ja — —Dap 



(2a) 
(2b) 



with 



Dap = D^Sap + (L»|| - D±_)uic 



(3) 



a diffusion tensor that incorporates the anisotropy of translational diffusion of elongated objects, with D\\ > and 
Dji the rate of rotational diffusion. The currents given in Eqs. (j2ap and \2h\ incorporate both diffusion and binary 
interactions. The latter are described by a mean field excluded volume potential. Vex, given by 



|U1 X U2|c(ri +S, U2,t) , 



(4) 



with m the mass of a rod, T the temperature, and s = r2 — ri the separation between the enters of mass of the two 
rods when they are at contact (Fig. ). If the thickness of the rods is negligible compared to their length, we can 



2 




FIG. 1: (color online) Long, thin rods of length i are parametrized by the position of their center of mass and a unit vector 
Ui denoting the rod's orientation in the plane. The position along the i-th rod measured from the center of mass is denoted by 
Si, with —ijl < Si < £/2. The overlap situation shown in the figure requires ri + uiSi ~ r2 + U2S2, so that s = siUi — S2U2 
The area excluded by rod 2 to rod 1 is the area of the dashed parallelogram spanned by the unit vectors ui and U2, given by 

£'\U1 X U2|. 



approximate s ~ siUi — S2U2, where —£/2 < Si < £/2 for i = 1,2, parametrizes the position along the i-th rod of 
length £ measured from its center of mass. The integral over the vector s spans the area excluded to rod 1 by a second 
rod with center of mass at r2, oriented in the direction U2. The excluded volume potential represents the second virial 
coefficient of the static structure factor and was first derived by Onsager [isj . 

In this paper we present the derivation of a modified Smoluchowski equation that describes the low density, over- 
damped dynamics of a collection of self-propelled hard rods. We consider long, thin hard rods moving on a substrate 
characterized by a friction constant C- Self propulsion is modeled by assuming that each rod moves along one direction 
of its long axis with a constant speed vq. In addition, the rods experience binary hard-core collisions. This is the sim- 
plest model for a "living nematic liquid crystal" , a terminology that has been used to describe the collective behavior 
of a variety of intrinsically sclf-propcUcd systems, from bacterial suspensions to monolayers of vibrated granular rods. 
Since the rods have a purely dynamical self-replenishing momentum, the statistical mechanics needs to be derived 
from the underlying trajectory dynamics. The details of the derivation are described in this paper. The outcome is 
a modified Smoluchowski equation for the configurational probability density of the form given in Eq. ([T]) , but where 
the translational and rotational currents acquire additional contributions due to self propulsion and take the form 

Ja^VQUiaC- DlpdipC- j-^{dipVe:r)c- ^^^^ , (5a) 



7^lc + -^(7^ly« 

Kb! 



Rl^lSP , (5b) 

2kBTa ^ ' ^ ' 
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where 

D^P = Da/3 + Ds UlaUip 

= D^5al3 + {D\\ + Ds - D±_)uiaUip , (6) 

with Ds = Wq/C, is the diffusion tensor. Self-propulsion modifies the familiar Smoluchowski equation for hard rods in 
several important ways. The first modification is the eonvective mass flux at the self-propulsion speed vq along the 
axis of the rod, described by the first term on the right hand side of Eq. (|5a)) . Secondly, self-propulsion enhances the 
longitudinal diffusion constant Dy of the rods, according to Z?|| — ?> Z?|| + Ds = D\\{\ + mv^/kBTa), as shown in eq. (|6]) 
This enhancement arises because self-propelled particles perform a persistent random walk, as recently pointed out by 
other authors [20| . Finally, the momentum exchanged by two rods upon collision is rendered highly anisotropic 
by self-propulsion. This yields the additional collisional contributions to the excluded volume interaction described 
by the last terms in Eqs. ([5a|) and (|5b[) . The precise form of these contributions can be found in Section IV. B. The 
Smoluchowski equation for self-propelled hard rods is the central result of this work. It has also been shown by us that 
the novel terms arising from self-propulsion have important consequences for the long-wavelength, long-time behavior 
of the system by introducing new terms in the coarse-grained equations for the dynamics of conserved quantities 
and broken symmetry variables. These hydrodynamie signatures have been reported in earlier work [l2| . Finally, an 
additional result of the work presented here is a purely dynamical derivation of the familiar Onsager excluded volume 
potential for equilibrium hard rods, given in Section IV. A. 

The layout of this paper is as follows. In Section II we analyze the trajectory of hard rods moving on a substrate in 
two dimensions. Using the fact that their trajectories are piece-wise differentiable, with singularities at the time of each 
collision, we derive an expression for the collision operator governing the momentum exchanged in a binary collision. 
In Section III we derive a formal hierarchy of Fokker-Planck equations governing the noise-averaged dynamics of a 
collections of self-propelled hard rods. In section IV we consider the limit of high friction with the substrate that 
yields a fast relaxation of the linear and angular momentum degrees of freedom, relative to that of the configurational 
degrees of freedom. This approximation, together with a low density closure of the Fokker-Planck hierarchy, allows 
us to derive the Smoluchowski equation. We conclude with a brief discussion. 



II. BINARY COLLISION OF HARD RODS 



Our model is a collection of self-propelled thin hard rods of length H. and mass m, confined to two dimensions. 
Although we will focus below on the limit of long, thin rods, to describe a binary collision we need to incorporate 
their finite thickness of the rods. We model each rod as a capped rectangle of uniform mass density, consisting of a 
rectangle of length i and thickness 2R ^ capped at the two short sides by semicircles of radius i?, as shown in Fig. 
^ . Each rod is described by the position r of its center of mass and a unit vector u = cos 0x -I- sin Oy directed along 
its long axis. The rods have head/tail, i.e., nematic symmetry. This symmetry is broken by self-propulsion that is 
implemented by assuming that a force F of constant magnitude and directed along the rod's long axis acts on the 
center of mass of each rod. The direction of the self-propulsion force will be referred to as the "head" of the rod and 
the unit vector u is chosen to point in the direction of the head, so that F ~ Fvl. 

The rods move on a passive medium that provides frictional damping to their motion. Any other physical or 
chemical process that may be present in the system is assumed to occur on a fast time scale, such that it can be 
modeled as an additive Markovian white noise. The dynamics of a self propelled rod is then described by Langevin 
equations for the center of mass velocity v = dtv and the angular velocity us = iuj = zdtO, where z is normal to the 
plane of motion. The equations of motion are given by 

mdtVa = -CapVp + FUa + TJa {t) , (7a) 

dtLo = -C'^OJ + m (t) , (7b) 

where (^ap = CW'^-a'^p + C-L [Sap — UaUp) is a translational friction tensor with (^y < reflecting the fact that frictional 
damping is smaller for motion along the long axis of the rod, Ci? is a rotational friction constant, and r] and 77/f are 
white noise terms with zero mean and correlations 

iVa it) VP (t')> = ^o.p (fl) S{t- t') , (8a) 

{VR{t)VR{t'))^^R5{t-t') . (8b) 

For simplicity, we assume that the noise amplitudes Aq,^ and have the same form as in equilibrium. 



Aa^(u) = 2kBTaUp{vL)/m , 

Afl = 2kBTaCR/I , 



(9a) 
(9b) 
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FIG. 2: (color online) Each self-propelled particle is modeled as a capped rectangle of uniform mass density, consisting of 
a rectangle of length i and thickness 2R capped at the two short sides by semicircles of radius R (the width of the rod is 
exaggerated for clarity). Self-propulsion is provided by a force of constant magnitude F directed along the "head" of the 
rod. The vector ^ locates points on the boundary of the rod relative to the center of mass. The unit vector is defined as 
= z X u, where z points out of the plane. 



with Ta an active temperature and / the component of the moment of inertia tensor I of the rod along the long axis, 

The last approximate equality in Eq. (|T0)) holds in the limits £ ^ 2R of long, thin rods. The active temperature Ta 
is generally different from the thermodynamic temperature of the system and is a measure of the noise amplitude. 

The rods interact vifith each other exclusively via hard-core interactions. The collisions are instantaneous and 
conserve energy and momentum of the colliding rods. To incorporate these interactions in the Langevin equations for 
the particles, we need to construct a collision operator that generates the instantaneous collision. We consider two rods 
and denote by t = the origin of time. The two rods travel freely with linear and angular velocities xi = (vi, Wi) and 
X2 = {'^'2,^2) until a time r(ri, when they come into contact, where = {r^, u^, v^, u)} is the phase point of each 
rod. Denoting by a prime the post-collisional velocities, x'^ = (v', u;,9i=i_2, the time dependence of the observables Xi 
can be written as 

X, (t) = x,Q (r (Fi, F2) -t)+ x[Q {t - t (Fi, F2)) . (11) 
The equation of motion for Xi is then 

9t.T,; = Ax, ,5(i-T(Fi,F2)) , (12) 

with Acci = Xi — x'^. The post-collisional velocities are easily calculated by imposing conservation of energy and of 
linear and angular momentum, = mv^ and = luji , with the result [isl . [l6| 

p- = Pj - LjiA , (13a) 
L', = L, - I, X k,,A , (13b) 

where kji is the unit normal at the point of contact of the two rods and is directed from rod j to rod i, and £,i is the 
vector from the center of mass of rod i to the point of contact of the two rods, as shown in Fig. |3l The magnitude of 
the momentum transfer A is 

A = "^^^21 ■ V12 ^^^^ 

1 + 2^ f^i X k2r) • • f^i X k2i) U2 X k2r) • • U2 X 



with V12 = V12 + tL'i X — (x'2 X ^2 the relative velocity of the two rods at the point of contact, and V12 = vi — V2. 
To render Eq. (jl2[) explicit we need an expression for the time r(Fi,F2). The condition of contact can be written as 
the requirement that ri(T)-|- ^i(t) = Y2{t) + (,2{t) for some value of ^1 and ^2 on the surface of the two rods, where 
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FIG. 3: (color online) A cap-to-side collision of two self-propelled hard rods (the width of the rod is exaggerated for clarity), k 
is a unit vector from rod 2 to rod 1 normal to the point of contact. Points on the side of the rods are identified by vectors ^i. 



Ti (r) = Vi + ViT and Ui (r) = cos {9i + ojir) x + sin {6i + utir) y. It is apparent from Fig. [3] that ri2 + ^1—^2 niust lie 
along k2i. The condition of contact can then be written as two scalar equations, given by 



(ri2 + $i-€2) 
(ri2 



k2i = |ri2+$i-$2l = 0, 
-^2) • X k2l) = , 



(15a) 
(15b) 



where all variables are evaluated at time r. Equations (|15ap and (jl5bp arc implicit equations for r(ri,r2). The first 
condition, Eq. (jlSap . imposes that two rods be in contact at any point along their surface. The second condition, 
Eq. (fTSb)) . det ermines the precise point on the surface of each rod. As an illustration we consider the collision shown 
in Fig. [3] when the cap of rod 2 comes into contact with a side of rod 1. In this case, ki2 = u^*" = z x ui and the 
surface of the rod 1 is parametrized by ^1 = siUi — Ru'ti with = z x Uj. The point of contact on rod 2 is simply 
^2 = |u2 + Rui- The contact conditions Eqs. (|15ap and (|15b|) become 



ri2 • - ^U2 • - 2i? = , 



ri2 ■ Ui - -U2 • Ui + Si 







(16a) 
(16b) 



Eq. (jl6ap requires the rods to be at contact at any point along the side of rod 1, while Eq. (|16bp determines the 
position of contact along the side of rod 1. 

To obtain an expression that can be used to eliminate r(ri,r2) from the equation of motion, Eq. (|12p we use the 
identity 



S (/ (x) - /o) 
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{x - Xo) 




9/1 






Ox la:o 





(17) 



Using Eq. (IT7|) . together with 



dt 



(ri2 + ^1-^2) ■ k2l|t=r = V12 (r) • k2 



(18) 
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we can rewrite the temporal 5 function in the equation of motion, Eq. (jl2|) as 



Sit^T (Fi, Ta)) = S (ri2 + ^-^2) ' ^21 



V12 • k2i 



(19) 



Finally, although Eqs. P5alll5b|) determine the contact condition, an actual collision results from the contact only if 
the pre-coUisional velocities at the point of contact are such that the two particles are moving towards each other, 
i.e., V12 • k2i < 0. Putting all of these results together, the equation of motion for the observable Xi can be written 
in the form 



dtx,=Til,2)x, , 
where T(l, 2) is a binary collision operator, given by 



(20) 



r(i,2) 



Ci,€2 "'k 



e -V12 • k2i 



V12 • k2i 



S (ri2 + ^1-^2) • k2i U (ri2 + $1-^2) ■ z X k2i 



(21) 

with 612 a substitution operator such that bi2Xi = x\. The integration in Eq. (PT|) ranges over all physical collision 
geometries. In the following, we will focus on the case of rods of large aspect ratio it ^ 2i?). In this limit, cap-on-cap 
collisions are rare relative to to cap-on-sidc collisions and will be neglected. For cap-on-side collisions, the J^^ /j^^^ 
can be given an explicit representation of the form 



'5 52 + 7: + '5 52 



(5 Fco„t(si,S2,k2l) 



+ (5 Si 



^ 

2 

(5 ( £21 - u:^ 



(5 k2i - U]^ + (5 k2i 



5 k21 + U.^ 



(22) 



where ■•■ = ds... and the last line simply defines a compact notation for the condition of contact. Equation (j2ip . 



with the expression given in Eq. (|22p for the range of integration, is the generator of collisional dynamics that will be 
used in the rest of this work. 

The above considerations are readily generalized to a system of N rods by considering only binary collisions since 
collisions among particles with hard core interactions are instantaneous and the probability of three or more particles 
being at contact at the same instant is of measure zero. In addition, the dynamics arising from the white noise due 
to the interaction with the substrate is described by a continuous generator and does not lead to any additional 
singularity. The derivation of the collision trajectory can then be carried out as before. The dynamics of a system of 
N self-propelled hard rods moving on a passive substrate in the A^-rod phase space F = {Fi,F2, ...,F7v} is controlled 
by a set of coupled Langevin equations for the linear velocities dtVi = and angular velocities u>i = x dtUi = zw^, 
given by 



mdtVi 



m^T -Cf^i 



Fn, 



(23a) 
(23b) 



where Ci is a friction tensor with components Ci,ap = ( C|| ~ C±)''^ia'Uii3 + C±^ai3, Vi ^-^id rj^ are Gaussian random forces 
with zero mean and correlations as defined in Eqs. (|8all8 bp . Noise sources associated with different values of the rod 
index i are uncorrelated. 



III. NON-EQUILIBRIUM STATISTICAL MECHANICS 



We are interested in the collective behavior of self-propelled rods in the limit of large friction with the substrate. 
In this limit one can consider a description that applies on time scales t 3> '^C|~^ that neglects the fast dynamics of 

the linear and angular velocities and only considers the dynamics of the coordinate degrees of freedom, {ri,Ui}^j^. 
The derivation of the overdamped dynamics must, however, be carried out carefully when dealing with the singular 
limit of hard particles, where the interactions, even though conservative, depend on the velocities of the particles. 
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With this goal in mind, we consider an observable of the system A{T), where F is an N-rod phase point. Using 
Ito calculus [l7[ , the stochastic equation of motion for the observable A can be derived from the equations of motion, 
Eq. (|23a|) and (|23b| . for the phase space variables, with the result 



OA (r,t) 

Ft 



LA (T,t) 



(24) 



where the operator L is given by 



L 



N 

E ■ ^ 

4=1 

m 



(25) 



where 7^^ = Ui x 5^. . The binary substitution operator bij contained in T{i,j) ^ T{Ti,Tj) replaces the velocities 
of the pair with their post-collisional values, according to bijA{vi,\-j,u}i,u}j) = A (v^, , o;^) , leaving the 

velocities of all other particles unchanged. Equation ([M)) describes the stochastic dynamics of any observable A(r) 
for given initial conditions in phase space. The quantity of interest here, however, is not the stochastic observable 
itself, but rather its ensemble averaged value for a given ensemble of initial conditions, pN (F), i.e.. 



Ait)) = /rfFp^(F)i(F,i) . 

/ e7is J 



(26) 



Equivalently, we can treat the phase space density as the dynamical quantity and write the ensemble average of an 
observable as 



A{t)) = f dTpN{r,t)A{r) 

I ens J 



(27) 



Equation (|27l) defines the dynamics of the phase space probability density, pj^{T,t). Taking the time derivative of 
both Eqs. ([25]) and (P7)) and defining an adjoint operator C, 



dV dtp{V,t)A[T) 



dV pN {T)dtA{T,t) 
dT PN {T)LA (F,0 



dV 



CpN{T,t) A{T) 



we obtain a Liouville-like equation describing the time evolution of of the phase space probability density, 

d 



dt 



+ £ Ujv (r, i) = , 



(28) 



(29) 



where 



L 



F ^ 
m 



N 

i=l 



(30) 



The single-particle part of C is identified by a simple integration by parts. To determine the binary collision operator 
T(i,j) one needs to explicitly construct the restituting collision, with the result 



r(i,2)p = 



ST 



Si,S2,k 



1 



Vi2-k 



e -Vi2-k p, 



(31) 
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where is the generator of restituting cohisions, i.e., b^j^A = A{xi,Xj), i.e., it replaees the post-collisional 

velocities of the pair with their pre-cohisional values, and 6 (Vcont ^^ij 52, , defined in Eq. (|22p . enforces the 
condition of contact. 

Finally, we average over the noise and define pjv = (p). The dynamical equation describing the evolution of the 
noise-averaged density is 

^ + c]pNir,t) = o, (32) 



where C is the Liouville Fokker Planck operator, given by 

i—l ^ ' 

-^Af^9„.„a.,, - Ia^OIj] - i ^ T(^,J) . (33) 

The formulation just described is exact and can be used, for instance, to evaluate time correlation functions for the 
system. To proceed, we will now restrict our attention to a low-density collections of self-propelled rods. In this case 
one can make systematic approximations to obtain the effective coarse-grained theory in the form of a Smoluchowski 
equation. 

IV. LOW DENSITY EFFECTIVE STATISTICAL MECHANICS 

In order to carry out systematic approximations in the low density limit, it is convenient to define a hierarchy of 
reduced distribution functions as 

fraiTi,...,Tra,t) = V"' j dT„,+l...drN PN(T,t) . (34) 

The Liouville Fokker-Planck equation, Eq. (j32p . can then be rewritten as an infinite hierarchy of equations for the 
reduced distribution functions. The m-th equation for /„ couples to fm+i- The resulting Fokker-Planck hierarchy is 
analogous to the BBGKY hierarchy for Hamiltonian systems and forms the starting point for constructing approximate 
theories to describe the system in various regimes of interest. 

At low density, we consider the first equation of the hierarchy for the one-particle distribution function /i (Fiji), 
given by 



dh{Ti,t) 



vhiT^t) ^ y"dr2 r(i,2)/2(ri,r2,0 , (35) 



dt 

where the one-particle operator 2? is given by 

= Vl•VrJl+a;l•7^l/l + -^II•VvJl--Ca0a,,„ (i;i^j/i)-C^a„, (c.i/i)--LAi^a,,„a„,^/i-iA^92^/i . (36) 

Equation (|35l) is a generalized Fokker-Planck equation that includes binary collisions. We are interested in the limit 
of large friction, where linear and angular velocities relax on fast time scales. Our goal is to construct an approximate 
closed equation for a local concentration field c (r, u, i), defined as 

c(ri,ui,i)=/ /i (ri,vi,ui,a;i,t) . (37) 



Specifically, we seek to derive a kinetic equation for self propelled particles that is analogous to the mean field 
Smoluchowski equation for thermal nematics. In the remainder of this section we present a systematic method for 
deriving such a closed kinetic equation for c in the limit of low density. 

We adopt the simplest phenomenological closure of the Fokker-Planck equation ([55]) of the form ji = 
fi (ri,<)/i (r2,t). With this closure, Eq. ([55]) becomes a Boltzmann- Fokker-Planck equation for the one-particle 
distribution function, 

j dV^T (1,2) h{T,,t)h{T2,t) (38) 
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To derive the Smoluchowski equation, in addition to the concentration field c(r, u,t) given in Eq. p7p . we introduce 
translational and rotational currents J"^ and J^ , defined as velocity moments of the one particle distribution function, 
with 



J^(ri,ui,i)= / vi /i (ri,ui, vi,u;i,i) , 

J vi 

J^'(ri,ui,t) = / u;i /i (ri,ui,vi,u;i,t) . 

J vi 



(39a) 
(39b) 



The dynamical equations for the concentration c and the translational and rotational currents are obtained by taking 
moments of Eq. ([38)) . with the result. 



+ — jJ - — UiaC+ dip {VaV 0)0 + 11113 {LUaVl3)c = ~Ia 

at m ^ m 



at 



where 



VlaVij3 fi{Ti,t) 



{uJaVfj) = - / L0iaVif3 fl{Ti,t) 



(ujaUJp) = i / wiQa;i/3 /i(ri,t) 
are the second moments of /i with respect to the translational and rotational velocities, and 

il = j dT2 [r(i,2)«ij/i(ri,t)/i(r2,t) 
/f = j dT2 [T(i,2)c^i„]/i(ri,i)/i(r2,0 



(40a) 
(40b) 
(40c) 

(41a) 
(41b) 
(41c) 

(42a) 
(42b) 



are the linear and angular momentum transfers, respectively, due to collisions with the other particles in the system. 

The equations for the translational and rotational currents contain frictional damping and relax on time scales of 
order m/C|| due to the interaction with the substrate. On time scales t 3> ™C|]~^ the fluxes can be approximated as 



lim J„ 



t>>m/C|| 



p 

-— Ul^C + dii3{VaVi3)c + 7lii3{u}aVp)c + 



^J[im_^ = -Cfl^ [dip(ujaVp)c + nip{ujaUJi3)c + 



(43a) 
(43b) 



In this regime we can fully describe the system's dynamics in terms of the one-particle configurational probability 
c(r, u, t). To obtain a closed equation, we must express the currents as functionals of the configurational probability. 



A. Thermal Hard Rods: Derivation of the Onsager Excluded Volume Interaction 

Before deriving the Smoluchwski equation for self-propelled hard rods, we show in this section how this method 
can be implemented to provide a derivation of the Smoluchowski equation for thermal hard rods, with the well-known 
mean field Onsager excluded volume interaction. We set the self-propulsion force F — and assume the noise is 
thermal, i.e. Tq = T. 

For large friction the relaxation of the linear and angular velocities of the rods is controlled primarily by the 
interaction of the rods with the substrate, rather than by interparticle collisions. We can then assume that the for 
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times t >> m/C|| the velocity distribution has relaxed to its noninteracting value. In the absence of collisions, the 
Fokker-Planck equation, Eq. ([55]) , for the one-particle distribution function takes the form 



Cj 



-^(Ja,? - I 1 ^ 



dvi^vip H dy.^dy^ ] fi 

m 



(44) 



where we have used the form of the noise amplitude in Eq. (|8all8bp to eliminate it in favor of the temperature T. 
The first term on the left hand side of Eq. ([^^ can be neglected for t >> m/C||. The solution of the non-interacting 
Fokker-Planck equation can then be written in a factorized form as 



with 



/i (ri,vi,ui,u;i,t) ^ c{ri,ui,t) fM{vi,uji) 



fM{vi,uji) = Aexp (~mvl/2kBT) cxp{^ lLj'f/2kBT) 



(45) 



(46) 



and A a normalization constant. In other words, we assume that on the time scales of interest the linear and angular 
velocity distributions have relaxed to their equilibrium forms. Using this expression, the velocity moments defined in 
Eqs. (|41aj|42b|) can be immediately calculated with the result. 



kr,T, 



-5al3 , 



m 

{VaUJf}) = , 
{uJaLOf}) = Sap 



(47a) 
(47b) 
(47c) 



Next, we need to evaluate the coUisional transfer contributions, /J and defined in Eqs. (|42a|) and (j42b[) . This 
requires calculating the mean linear and angular velocity transferred in a collision, (T (1, 2) via) and (T (1, 2) ujia) 
where 



{X) 



M 



X /a/ ivi,UJi) /a/ (w2, W2) 



(48) 



Using the explicit form of the momentum transfer in a binary hard rod collision given in Eq. (|14p. we find 



(T(l,2)vi„) 



A/ 



Vi2-kl e(-Vi2-k 

SlS2fe \ ^ 

1 

1 + f (|ixk)' + f (6xk 



M 



(j'cont (^Sl,S2,k^) fco 



(49) 



The 6 function in Eq. (|49| selects only those values of the pre-collision velocities that will actually result in a collision. 
In the mean field limit of interest here, we assume that on average half the particles in the flux will have initial velocities 
that will yield a collision and approximate the velocity average in Eq. (j49p as 



Vi2-k e -V 



M 



M 



Substituting Eq. ^ in Eq. (gH) we find 

(r(i,2)vi„)^. 

Similarly, it is easy to show that 

(T(l,2)wi„)^, = - 



knT 



2kBT ksT 



xk 



.'\2 kuT 



6xk 



<5 r, 



Sl, S2,k) ) ka . 



ksT 



S1S2 -^k 



i^cont (si,S2,k)) (s X k^ 



(50) 



(51) 



(52) 
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where s is defined in Fig. [T] In the thin rod hmit, i.e., i? — >■ 0, the contact delta function in Eqs. ([5T|) and ((52)) is non 
zero on the perimeter of a parallelogram centered at the position ri rod 1 with sides of length £ directed along Ui and 
U2 (Fig. [1]). The area of such a parallelogram can be written as 

A||g^<,„(ri, ui) £2 / e (0+ - |ri2 + siui - S2U2I) • (53) 

J S1,S2 

It is easy to verify that 

VriA||g,.Q,„ = j 5 (^Fcont (^Sl,S2,k21^^ £21 , (54a) 

^l^llgram = / / (5 (Fcont (si , S2 , k) ) (s X kai) ■ (54b) 

Finally, combining all these results we obtain the familiar Smoluchowski equation for hard given in Eqs. ([l}, (|2bp and 
(|2ap with Dafj = Cafi^BT/m and Df; = ksT / {C^nl). The excluded volume interaction can be written in the familiar 

Onsager form by noting that the area of the coUisional parallelogram is simply given by A^^^^"-™" = x U2I as 

illustrated in Fig. ([1]). This gives Eq. where the integral over s spans the area of the parallelogram shown in 
Fig. ([T]) for fixed orientation of rod 2, while the integral over U2 averages over all possible orientations of the second 
rod. The Onsager excluded volume interaction has a purely entropic interpretation as the free energy cost for a rod 
to occupy the area excluded by another rod. Here the same result has been derived from dynamical considerations 
and has the alternate interpretation of mean-field momentum transfer in a binary collision of two rods, each carrying 
an average momentum of magnitude y/rnksT. The oricntational correlations arise then from the anisotropy of the 
collision frequency due to the fact that head to head collisions are of measure zero with respect to head to side 
collisions. 

In summary, we have shown in this subsection that for thermal (non self-propelled) hard rods the method described 
in this paper can be used to derive the familiar Smoluchowski equation with the Onsager expression for the excluded 
volume interaction. In the next section we apply the same procedure to self-propelled hard rods and show that 
self-propulsion modifies the Smoluchowski equation in several ways. 



B. Self-propelled hard rods 



As for the case of thermal hard rods, the derivation of the Smoluchowski equation for self propelled rods consists 
of two steps: (i) we identify a stationary normal solution of the noninteracting Boltzmann-Fokker Plank equation, 
and (ii) we use this functional form to close the equations for the translational and rotational fiuxes in a quasi-static 
approximation. 

In the absence of interactions, the Fokkcr-Planck equation for self-propelled rods is given by 



Cii 



{dt + vi • Vri -(- u;i • 7^l) /i -I- wqUi-Vvi 



Sap 



1 



Cii 



m ' 



knT. 



(55) 



with vq = the self propulsion velocity. On time scales large compared to m/C||, we neglect the first term on the 
left hand side of Eq. . The stationary normal solution then has the form 



fl (ri,ui,vi,a;i,t) = c(ri,Ui,t)/s (vi,a;i|ui) , 



with 



fs^A exp 



-m (vi - vqUi) 



1 



-Ilu7 



1 ' 



2fcsT ' ' 2kBT 

and A a normalization constant. Inserting this stationary normal solution in Eqs. (|43l) . we obtain 



JI =VoUlaC- DsUiaUl-Vr^C- Dafjdr^^C-mC^p (^j dT2dVL2 (^(1,2 



) vi/3)oc(r2,U2,i) c(ri,ui,i) 



= -DbUic -^{^J dr2du2 {T (1, 2) c (r2, U2, t)^ c (ri, ui, t) , 



(56) 
(57) 

(58a) 
(58b) 
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where the angular brackets {■.■)s denote the average over the hnear and angular velocities with weigh fs- The diffusion 

2 

coefficients Dap and are as given in the previous subsection, but with T ~Ta, and Ds ~ 

The next step is the evaluation of the average linear and angular momentum transfer in a collision, 



(r(l,2)vi)5= / / T(l,2)vi /s(vi,u;i|ui)/s(v2,a;2|u2) , 

J Vi ,UJ\ J V2 ,C*^2 

(r(l,2)u;i)5- /" I T(l,2)a;i /5(vi,u;i|ui)/5(v2,u;2|<i2) . 



(59a) 
(59b) 



When deriving the Onsager excluded volume interaction in the previous subsection, we neglected dynamical velocity 
correlations by approximating (^(Vi2 • k)^8(— V12 • k)^ ~' i ^(^12 • k)-^^ This approximation is not valid for self- 
propelled rods as the dynamical velocity correlations induced by the self propulsion velocity, which is directed along u, 
are also orientational correlations. On the other hand, an approximation is required to make the problem tractable. 

The Smoluchowski equation is a mean field model and only describes the average momentum transfer in a collision. 
We then assume 



{T (1, 2) vi)^ ^ (T (1, 2) vi)^,^ + (T (1,2) v^^p , 
(T(l,2)a;i)^ ^ {T (1,2) u^) + {T (1,2) lj,) , 



(60a) 
(60b) 



where {X) denotes velocity averages with the Maxwellian distribution given in Eq. (pS]) at temperature T = Ta and 
(X) gp denotes velocity averages in a regime where all particles are moving at velocity vqu and the one-rod velocity 
distribution is given by fsp ^ S {v — vqu) 5 {u:) . The first terms on the right hand side of Eqs. are then evaluated 
as in the previous section by neglecting the dynamical velocity correlations and yield again the Onsager excluded 
volume potential given in Eq. The second terms on the right hand side of Eqs. (|5n)) are evaluated by neglecting 

r - - 1 

static orientational correlations arising from the term 1 + (m/2/)(^i x k)'^-|-(m/2/)(^2 x k)^ in Eq. as these 
correlations are already incorporated in the noise. It can be shown that this approximation becomes exact if we 
replace the rod by a string of beads in contact with each other and calculate the momentum transfer between the two 
specific beads that participate in the collision. With these approximations, the average momentum transfer is given 
by 



1 



2 tSP 



with 



rSP 



{T (1, 2) vi^)^ c (2) j c (1) ~ - (Vr,„14.) c + vll[ 
^(r(l,2)a;i)5c(2)) c(l) ^ i7^lK.c+ , 



sin^ (6*1 - O2) [e (sin (6*1 - B^)) - 6 (- sin {Oi - 62))] 



(61a) 
(61b) 



Ui c ( ri + sui - -U2,U2,t I + U2 c I ri + SU2 - -Ui,U2,i 



(62) 



fSP 



S C\Vi 



h - 62) [0 (sin (01 - 62)) - e (- sin (01 - 92))] 



SUi - -U2,U2,t 



cos {9i - ^2) c ri + SU2 - -ui, U2, t 



(63) 



Finally, when these results are substituted into Eq. (|58ap and (|58bp . we obtain the modified Smoluchowski equation 
given in Eq. ([5]). 

There are three important modifications of the Smoluchowski equation for self-propelled particles when compared 
to its equilibrium counterpart for thermal particles. 

1. The translational current in Eq. ()5ap contains a convective term vqvlic that arises because self-propelled 
particles move in the direction of their long axis. This is a signature of the polar nature of the microdynamics 
of the system. 
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2. Self propulsion yields an additional longitudinal diffusion current DsuiaUipdipc, with D5 ~ Vq. This can be 
understood as follows. A Brownian particle subject to a frictional damping C takes Brownian steps of mean 
length A ~ \/mkBT/C., with ^JkBT/m the thermal speed of the particle and m/C the frictional time scale over 
which the velocity relaxes to zero. This yields the simple estimate D ^ A^/(m/C) ^ ksT/C, for the diffusion 
coefficient. A Brownian rod experiences anisotropic friction C|| and yielding mean steps A|| ~ y/mkBT / C,\\ 
and Ax ^ \JmkBT / C,i_ in the directions longitudinal and transverse to its long axis, respectively. This gives 
anisotropic diffusion constants D\\ ^ kBTC,\\ and D± ^ kBT/(±. When the rod is self propelled it performs a 
persistent random walk along its long axis controlled by the competition between ballistic motion at speed vq and 
rotational diffusion at rate ^ _D||/£^. As a result, the mean square velocity is approximately (v^) ^ ■^^^ + ^0- 

For small values of vo this gives {v^)'^/'^ ^ 
to the longitudinal diffusion, corresponding to the second term of Eq. 

3. Both the translational and rotational fluxes in Eq. ([5]) contain additional active contributions arising from 
the momentum transfer associated with self propulsion. Collisions among particle induce dynamical velocity 
correlations. These are negligible for a mean-field description of overdamped thermal hard rods that only seeks 
to capture the dynamics of the translational and orientational degrees of freedom. When the rods are self- 
propelled with a physical velocity directed along their long axis, collisions induce both velocity and orientational 
correlations since the two are intimately coupled. These additional collision-induced orientational correlations 
have been shown to affect the physics of the system even on hydrodynamic scales (l2j . 

Our work demonstrates that orientational fluctuations have a profound effect on mass transport in self-propelled 
particle systems. Furthermore this observation is not limited to the specific hard rod model considered here, but holds 
generically for all collections of self-propelled units. Our result relies solely on the presence of short-range excluded 
volume interactions that are present in all physical systems, but not on the specific implementation of such interaction 
in the hard particle limit. 

V. DISCUSSION 

In this paper we have analyzed the microscopic dynamics and statistical mechanics of a collection of self-propelled 
particles modeled as long thin polar rods that move along one direction of their long axis. The formalism developed 
here is general and of wide applicability. It can be used to study fiuctuations and response functions in fluids of self- 
propelled particles by building on techniques developed for traditional fluids. In addition, the formalism can readily 
be generalized to particles of arbitrary shape and to higher dimensions, making it also relevant for applications to 
granular fluids. 

As a particular application of the general formalism, we have derived the Smoluchowski equation that governs the 
dynamics of the one particle configuration probability density in the overdamped regime, for both thermal and self 
propelled particles. For thermal rods, this provides a purely dynamical derivation of the familiar Onsager excluded 
volume interaction and is useful for identifying the limitations of this widely used effective interaction in capturing the 
physics of out-of-equilibrium systems. For self propelled hard rods, the modified Smoluchowski equation presented 
here is the first tractable theory of the dynamics of self-propelled particles that captures the physics of contact 
interactions and their modifications due to the presence of self propulsion. In a separate work we have used this 
Smoluchowski equation as the starting point for deriving a coarse grained (hydrodynamic) description of the system 
in terms of conserved quantities and broken symmetry variables [l3 . [isj . Self propulsion has profound effects on the 
system on hydrodynamic scales, including enhanced orientational order arising from the orientational correlations 
induced dynamically via collisions, instabilities of the ordered phases, and the existence of propagating sound-like 
density waves in this otherwise overdamped system. 

Having derived the Smoluchowski equation from first principles, we can also identify its scope and limitations. Our 
theory captures the fluctuations in velocity induced by orientational fluctuations and the associated modifications 
to the mass fiux that characterize this inherently nonequilibrium system. These effects yield the convective and 
diffusive terms in Eq. ([5]). The theory also captures some of the orientational correlations induced by collisional 
interactions through the additional momentum transfer contributions to the fluxes. These correlations are important 
as evidenced by the enhanced ordering identifled in [l^ and observed in numerical simulations of motility assays 
0. The derivation is, however, based on a low-density kinetic theory that neglects two-particle velocity correlations. 
While this is a reasonable approximation for overdamped thermal particles that have an underlying equilibrium state, 
in the case of self propelled particles, these correlations will generate additional orientational correlations that are 
neglected in the present approximation. The Fokker-Planck hierarchy derived here can serve as the starting point for 
analyzing the effect of these two-particle correlations. This is left for future work. 



2kBT 



This yields an additional contribution of order 
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In summary, we have constructed the non-equihbrium statistical mechanics of a system of self propelled particles 
and derived a modified Smoluchowski equation for the system. We have discussed the content of the resulting theory 
and identified its scope, limitations and potential for future applications. 
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